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Abstract. We sketch here an approach to the computation of genaralised parton distributions 
(GPDs), based upon a rainbow-ladder (RL) truncation of QCD’s Dyson-Schwinger equations 
and exemplihed via the pion’s valence dressed-quark GPD, Our analysis focuses 

on the case of zero skewness, ^ = 0, and underlines that the impulse-approximation used 
hitherto to define the pion’s valence dressed-quark GPD is generally invalid owing to omission 
of contributions from the gluons which bind dressed-quarks into the pion. A simple correction 
enables us to identify a practicable improvement to the approximation for H^{x, 0, t), expressed 
as the Radon transform of a single amplitude. Therewith we obtain results for (x, 0, t) and the 
associated impact-parameter dependent distribution, q^ix, |6x|), which provide a qualitatively 
sound picture of the pion’s dressed-quark structure at an hadronic scale. 


1. Introduction 

The current contribution is aimed at reporting on the results, recently presented in Refs. mm, 
concerning the computation of the pion’s valence dressed-quark GPD, based upon a rainbow- 
ladder (RL) truncation of QCD’s Dyson-Schwinger equations (DSE), and its connection to the 
parton distribution function (PDF) in the appropriate forward limit. 

After more than forty years of studying PDFs, since quarks were discovered in a series of 
deep inelastic scattering (DIS) experiments at the Stanford Linear Accelerator Center [31 [HIS], 
a good deal is known about hadron light-front structure. Notwithstanding that, much more 
needs to be understood, particularly, e.g., in connection with the distribution of helicity 
iiE]. Moreover, PDFs only describe hadron light-front structure incompletely because inclusive 
DIS measurements do not yield information about the distribution of partons in the plane 
perpendicular to the bound-state’s total momentum; i.e., within the light front. Such information 
is expressed in generalised parton distributions (GPDs) (HHH [131 [H] , which are accessible via 
deeply virtual Compton scattering on a target hadron, T; viz., 'y*{q)T{p) —)• 'y*{q')T{p'), so 
long as at least one of the photons [7*((?), 7*(90] possesses large virtuality, and in the analogous 
process of deeply virtual meson production: 'y*{q)T{p) —)• M{q')T{p'). Importantly [see Sect. 2], 
GPDs connect PDFs with hadron form factors because any PDF may be recovered as a forward 
limit of the relevant GPD and any hadron elastic form factors can be expressed via a GPD-based 


sum rule. The potential that GPDs hold for providing manifold insights into hadron structure 
has led to intense experimental and theoretical activity naiiaiiaiiaiis]. 

Most of the constraints that apply to GPDs are fulhlled when the GPD is written as a double 
distribution mini US], which is equivalent to expressing the GPD as a Radon transform m- 

H{x,^,t)= f dadp 5{x — a — ^f3) [F{a, I3,t) + ^G{a, /3,t)] , ( 1 ) 

./|a|+|/3|<l 

where the variables x, t are defined following Eq. Q and, at leading-twist, F, G have 
operator definitions analogous to the GPD itself. In order to obtain insights into the nature 
of hadron GPDs, it has been common to model the Radon amplitudes, E, G, following 
Refs. [5U]. This approach has achieved some phenomenological success ffSlEI]; but more flexible 
parametrisations enable a better fit to data [22] . Such fits played a valuable role in establishing 
the GPD framework; but if one wishes to use measured GPDs as a means by which to validate our 
basic perception of strong interactions in the Standard Model, then data fitting is inadequate. 
Instead, it is necessary to compute GPDs using a framework that possesses a direct connection 
with QGD. This observation is highlighted by experience drawn from the simpler case of the 
pion’s valence-quark PDF [T]. Herein, we therefore adopt a different approach, sketching a 
procedure for the computation of hadron GPDs based on the example provided by the pion’s 
valence-quark PDF. 

2. General features of pion GPDs 

From a quark model perspective, in the isospin symmetric limit, the pion is a quantum 
mechanical bound-state of two equal-mass constituents and it is therefore the simplest hadronic 
bound-state. That is a misapprehension, however. Owing to the connection between pion 
properties and dynamical chiral symmetry breaking (DGSB); i.e., its dichotomous nature as a 
Goldstone mode and relativistic bound-state [23l[2l]) a veracious description of the pion is only 
possible within a framework that faithfully expresses symmetries and their breaking patterns. 
The Dyson-Schwinger equations (DSEs) fulhll this requirement [25112611^ and hence we employ 
that framework to compute pion properties herein. 

Notwithstanding the complex nature of the pion bound-state, it is still a J = 0 system and 
hence for a vector probe there is only one GPD associated with a quark q in the pion (vr^, tt^). 
It is dehned by the matrix element 

/ d^7 

— 6in • (vr(P^.)|g- {-z/2) n-^q {z/2) |7r(P_)), (2) 

where: k, n are light-like four-vectors, satisfying = 0 = n^, k ■ n = 1] z±^ represents that two- 
component part of z annihilated by both k, n; and P± = P± A/2. In Eq. ([2|), ^ = —n-A/[2n-P] 
is the “skewness”, t = — is the momentum transfer, and P^ = — P-A = 0. The GPD 

also depends on the resolving scale, </• Within the domain upon which perturbation theory is 
valid, evolution to another scale C' is described by the ERBL equations [28l[29] for \x\ < ^ and 
the DGLAP equations [Ml EH |32l EH] for |x| > where ^ > 0. 

In order to produce quantities that are gauge invariant for all values of z, Eq. ([2]) should 
contain a Wilson line, Wl—.2;/2,z/2], between the quark fields. Notably, for any light-front 
trajectory, >V[—z/2,z/2] = 1 in lightcone gauge: n ■ A = 0, and hence the Wilson line does 
not contribute in this case. On the other hand, light-cone gauge is seldom practicable in 
either model calculations or quantitative nonperturbative analyses in continuum QGD. Indeed, 
herein, as typical of nonperturbative DSE studies, we employ Landau gauge because, inter 
alia |34[ 135] : it is a hxed point of the renormalisation group; and a covariant gauge, which 
is readily implemented in numerical simulations of lattice-QGD. It is therefore significant that 


VV[—z/2, zl‘2\ is not quantitatively important in the calculation of the leading-twist contributions 
to numerous matrix elements |36j . 

It is worth recapitulating here upon some general properties of GPDs. Most generally, 
Poincare covariance entails that GPDs are only nonzero on x € (—1,1). Moreover, owing 
to time-reversal invariance, Kinematically, the skewness is bounded: 

^ € [—1,1], but ^ E [0,1] for all known processes that provide empirical access to GPDs. 

Focusing on the pion, a charge conjugation mapping between charged states entails 

and consequently, in the isospin symmetric limit: 

= (3) 


It follows that the isospin projections: 

H\x,i,t) :=F“+(x,C,t) + (-l)"^^+(x,?,t), 1 = 0,1, (4) 

have well-defined symmetry properties under x -fE —x; viz., is odd and is even. 

Returning to the definition in Eq. ([2]), it is plain that if one considers the forward limit: ^ = 0, 
t = 0, then X is Bjorken-x and the GPD reduces to a PDF; viz.. 




q^{x), X > 0 
—x), X < 0 


(5) 


Moreover, irrespective of the value of the electromagnetic pion form factor may be computed 
as 


F,+ (A2) = + (6) 

=: (A2) + (A^) = (A^) , (7) 

where are the quark electric charges in units of the positron charge and we have used Eq. ([3]) 
to show E^+(A^) = —E“+(A^). Additional information may be found elsewhere |37j . 


3. Heuristic Example 

Imagine a bound-state of two scalar particles with effective mass a and suppose that the 
interaction between them is such that it produces a light-front wave function of the form 
(x = 1 — x): 





'XX 


XX 


-6{x)6{x) 


( 8 ) 


(A merit of considering a bound-state of scalar constituents is that in describing the wave 
function of the composite system one avoids the complication of Melosh rotations, which arise 
in treating spin states in light-front quantum mechanics [38] •) If the skewness is zero, in which 
case the momentum transfer is purely light-front transverse, then the GPD for this system can 
be written as a wave function overlap |39llini[l2l[l3|: 


A 


(fk± tpix, k_i_ + (1 — x)A_i_) -(/^(x, k_i_). 


(9) 


This entails 


{77,(x,0,Ai) >0:xG [-l,l],Al >0}. 


(10) 





Owing to the simplicity of the starting point, Eqs. ([ 8 ]) and (l9|) allow one to obtain an algebraic 
expression for the GPD; viz., with ^ = A^(l — x)/4x(T^, then 


Hcr{x, 0, —A]_) = 30(1 — xyx"‘C{z) 6{x)6{x ), 


In 


C{z) = 


z^ + (2^ + 1) yz^ ^ \ + 32; j / ( + 4 — z 


+ 4 


( 11 ) 

( 12 ) 


Some further analysis reveals that C{z) decreases monotonically away from its maximum value 
C{z = 0) = 1. Consequently, Hfj{x,0,0) = 30(1 — x)^x^0(x)0(x), which is an excellent 
approximation to the pion’s valence dressed-quark PDF [T]; and whilst the maximum of 
Hfj{x,0, —A'^) lies at X = 1/2 for A^ = 0 , this peak shifts to x = 1 with increasing A^, 
consistent with an expectation that for A^ ;§> the interaction probability is largest when the 
probe and hadron are collinear | 12 ] . 

The Hankel transform: 

qaix,\by) = I ^|Ax|Jo(|6x||Ax|)^(^,0,-Ai), (13) 

dehnes the system’s impact-parameter-dependent (IPD) GPD [39]. It is a density that describes 
the probability of hnding a parton within the light-front at a transverse position b± from the 
hadron’s centre of transverse momentum (CoTM). Since i^(x,0, — A^) is a positive-dehnite, 
monotonically decreasing function of A^ for each x, the global maximum of ( 7 o-(x, | 6 x|) is located 
at | 6 _l| =0 and q^ix, | 6 _l|) is a monotonically decreasing, positive-definite function away from 
that maximum. 

The value of x at which the global maximum in q„(x,\b±\) occurs is determined by the 
system’s dynamics. Considering the hadron’s valence dressed-parton structure, one extreme 
is achieved if C{z) is independent of x: the maximum of qo-{x,\b±\) is then located at 
(x = 1/2, | 6 _l| = 0). In realistic cases, the necessary (x,A^) correlations in C{z) and this 
function’s general properties act to shift the maximum to x > 1/2. Using Eqs. ([ 8 |) and ([9]), 
qa{x,\b±\) peaks at (x = 0.72, | 6 _l| = 0). One may also consider the path followed by the 
maximum as one increases \b±\ away from zero. To that end, observe from Eq. (11311 that for 
|^±| S> l/cr the x-dependence of qa-{x, | 6 _l|) is dominated by H{x, 0 , —A^ ~ 0 ), which peaks at 
X = 1/2. The nature of C{z) then entails that the peak in the valence dressed-quark IPD GPD 
drifts monotonically toward x = 1/2 as cr\b±\ —)• 0 . 


4. Pion’s valence dressed-quark GPD 

In order to compute an approximation to the valence-quark piece of the GPD expressed in 
Eq. ([2]) we adapt the method used successfully elsewhere to compute the pion’s valence-quark 
distribution function [T] and elastic form factor |2]- Consider, therefore, 

2H:ix,tt) = N,tr [ (14) 

JM 

where := J is a translationally invariant regularisation of the integral; {£) := 

5{n ■ £ — xn ■ P); the trace is over spinor indices; 77 G [ 0 , 1 ], f/ = 1 — 77 ; £^ = fji^ + r}£p, 
_|_ fj£p^ £± = £± A/2, £p = £ — P. (N.B. Owing to Poincare covariance, no observable 
can legitimately depend on 77 ; i.e., the definition of the relative momentum.) So long as each of 
the dressed-quark propagators, S{£), on the right-hand-side (rhs) of Eq. (fTT|) is computed using 
the rainbow truncation of QCD’s gap equation, and both the pion Bethe-Salpeter amplitudes, 










r^(£; P), and the dressed-quark-photon vertex, are calculated in the associated ladder 

truncation of the relevant Bethe-Salpeter equations then t), thus computed and inserted 

in Eq. Q, provides the leading-order contribution to the pion’s electromagnetic form factor in 
the most widely used, symmetry preserving truncation of QCD’s DSEs: the rainbow-ladder (RL) 
truncation [421 143] , whose strengths and limitations are detailed elsewhere [251 ESI EZ]. 

Given its connection with a reliable scheme for computing it was long thought [441145j 

that Eq. (I14p would also be an adequate starting point for computation of the pion’s valence- 
quark PDF, qy{x). However, as explained in Ref. [T], that is not true. 



Figure 1. Amplitude-(l) for virtual-photon-pion Compton scattering in RL truncation is 
obtained from the sum (A) -|- [B) — (G). The “dots” in (A) and (B) indicate summation of 
infinitely many ladder-like rungs. The other two amplitudes are obtained as follows: (2) - 
switch vertices to which q and q' are attached; and (3) - switch vertex insertions associated 
with q' and P'. The lines and vertices mean the following: dashed line - pion; undulating line 
- photon; spring - interaction-gluon in the RL kernels; solid line - dressed-quark; open-cross 
cirele - dressed-quark-photon vertex; filled cirele - pion Bethe-Salpeter amplitude. Each of the 
last three is computed in RL truncation. 

Diagram-(A) in Fig.[T] represents the textbook handbag contribution to virtual Compton 
scattering and, in the forward and Bjorken limits, yields the following piece of the structure 
function: 

q\{x) = Ncti f 6f^{kn)n--/'H-,r{P,k) , (15) 

Jdk 

which is typically the only contribution retained in computations of the pion’s quark distribution 
function; e.g., Refs. HSlllZlllHl Ha El is]. In RL truncation, fc) is an infinite sum of 

ladder-like rungs, as illustrated in Fig.[Tl so that one may write |45| 

q\{x) = Ncti [ iTT,{kr^,-P) S{kn)T'^{k;x) S{kn)iT^{kfj,P) S{kf^), (16) 

Jdk 

wherein P^(k]x) is a generalisation of the quark-photon vertex, describing a dressed-quark 
scattering from a zero momentum photon. It satisfies a RL Bethe-Salpeter equation with 
inhomogeneity in ■ ^ 5f^^{krf) [IS]- Equation (fTUl) ensures fgdxq^(x) = 1 because f dxr^(£;x) 
gives the Ward-identity vertex and the Bethe-Salpeter amplitude is canonically normalised. 
The minimal Ansatz sufficient to preserve these qualities is T^{k;x) = 6^{krj) dk^S~^{krj), in 
which case one has 

q\ix) = Wtr [ 6f,{k^)r^{k^, -P)dk^S{k^)K{k^, P) S{k ^). (17) 

Jdk 

As discussed in Ref. |2], the forward limit for the GPD defined by Eo. (fT41) . L7^(x,0,0), 
corresponds exactly to the pion PDF contribution given by (fTTll . q\{x). However, as proven in 














Ref. [T] , that impulse approximation is incomplete because it omits a fraction of the contributions 
from gluons which bind dressed-quarks into the pion. These contributions appear depicted in 
the diagram of Fig.[2J As has been lengthly explained in Ref. [1], the contribution depicted 



Figure 2. Forward limit of the combination 
{B) — (C) in Fig.[T] The figure also 
exposes the internal structure of the pion’s 
Bethe-Salpeter amplitude obtained in RL- 
truncation. In the Bjorken limit, the 
initial/final-state interactions involve very 
soft gluons and hence, in combination 
with the ladder resummation, produce a 
contribution that is of the same order as 
Diagram-(A) in Fig.[TJ 



Figure 3. Imaginary part of the left panel 
in the Bjorken limit: the vertex insertion 
can appear between any two interaction 
lines. The compound vertex on the right 
is readily simplified using the RL Bethe- 
Salpeter equation. 


in Fig.O which has hitherto been overlooked, is of the same order as that from Diagram-(A), 
leading to the result of (1171) . Consider now the imaginary part of Fig. [2] in the Bjorken limit, 
which produces the leading contribution illustrated in the Fig.O the vertex insertion can appear 
between any pair of interaction lines. Using the recursive structure of the ladder Bethe-Salpeter 
kernel and the Ward identity, which entails that inserting a zero-momentum vector-probe into 
a propagator line is equivalent to differentiation of the propagator, then the compound vertex 
on the right side of Fig. [3] is readily seen to correspond to differentiation of the Bethe-Salpeter 
amplitude itself with respect to One thus arrives at the following contribution from (B) — [C) 
to the pion’s quark distribution function; 

qlcix) = N,tl [ 6^ikr,)dk,K{k^, -P)S{kr,)r^{k^, P) S{kf^) . (18) 

Jdk 

This expression is nonzero in general and only vanishes when the pion’s Bethe-Salpeter amplitude 
is independent of relative momentum, so that one is left with 

Qvix) = q\{x) + qlcix) , (19) 

for the total quark distribution in RL truncation. 

Therefore, since Eq. (fTTl) is incomplete for qy{x) it is necessarily also inadequate for 
computation of H^{x, t). Importantly, we have found that the flaw is expressed more forcefully 

as t := grows: one can obtain H^{x, 0,—t) < 0, which is physically impossible, as explained 
in connection with Eqs. dSD, m- 







The amendment to Eq. ()14H is to be guided by the precise knowledge of the form for the 
correction q^(^{x), as the former should be non-forward latter’s extension in the case of the pion 
valence-quark GPD. However, it is still being sought and below, for illustrative purposes, we 
content ourselves with the following ansatz to aument Eq. (fTTll = n • cI^r] 


^ Jm 

xS{l.) + -P+)5(£p)d!ir^(£^; P-)S{i.)] . 


( 20 ) 


This Ansatz extends the handbag diagram correction identified in Ref. [T] to t > 0; and, 

in connection with the valence dressed-quark GPD, it can be expressed via a Radon amplitude 
F^{a,l3,t) which preserves the good features of the kindred amplitude produced by Eq. (fTTl) . 

We will now focus primarily on = 0, be guided by Eq. (I14p . and mention and ameliorate its 
failings where appropriate, drawing on the insights gained from the example in Sect. 3. Notably, 
the defects of Eq. (|14l) are typically overlooked in extant continuum computations of the pion’s 
GPD [501 [5T1 [52l [53l[5l] ): Refs. |5T1[53] deliver a form for q{x, [ftpl) that is not positive dehnite. 

In order to gain novel insights into pion structure, we use the following algebraic forms for 
the dressed-quark and pion elements in Eq. (fTT|l + M^)] [55] : 


s{e) 

Pu{z) 


ru,T^{il-±P) 


H7-£ + M]Am(^'), 

1 r(u + 3/2) _ 2 

0F r(z/ + i) ^ ^ ’ 

W5 j ^dz p^iz) , 


(21a) 

(21b) 

(21c) 


where M is a dressed-quark mass-scale; Am(-^^) = M^Am(^^); ^z± = (^ + {zP 1)^/2 and we 
work in the chiral limit (P^ = 0 = m, where m is the current-quark mass); and is the Bethe- 
Salpeter amplitude’s canonical normalisation constant. Owing to the simplicity of Eqs. (1211) . one 
can reasonably employ r^(^+,.^_) = 7 ^Pr(t = A^), where Prit) is the vertex dressing function 
described in Eqs. (24)-(28) of Ref. [55] . 

Eirst, working with the input specified in connection with Eqs. (|2ip . one can obtain the two 
contributions for the valence-quark parton distribution, qyix), from Eqs. (I17ll8p . 

72 

Qa{^) = —ra:^(x[—2(x — 4)x — 15]-|-30) ln(x)-|-(2x^-|-3) 

25 

X (x — 1)“^ ln(l — x) -b x[x(x[2x — 5] — 15) — 3](x — 1)] , (22a) 

Qsci^) — ^ [ic^(2x([x — 3]x-b 5) — 15) ln(x) — (2x^-b 4x-b 9) 

Zo 

X (x — 1)^ ln(l — x) — x(2x — l)([x — l]x — 9)(x — 1)] ; (22b) 

and display them in Eig. [H Only after summing the two contributions, q\{x) and q^Q{x), 
considering hence the leading contributions from all the diagrams in Eig.[Tl the valence-quark 
parton distribution is symmetric under the exchange x —> 1 — x: viz. the dressed-quark and 
-antiquark each carry half the pion’s momentum. 

We computed then, with the same input specihed in connection with Eqs. ()2ip . the triangle 
diagram result from Eq. ()14p and the correction given by Eq. ()20p . both to obtain all the leading 
contributions to the pion’s valence dressed-quark GPD, Hi^{x,^,t). As detailed elsewhere [37] . 
that task for triangle diagram was completed by deriving an expression for the Mellin moments 
of H^{x,^,t) from Eq. (fTTl) and introducing five Eeynman parameters {x,y,u,v,w), dehned in 







Figure 4. Pion dressed-quark distribution 
function in rainbow-ladder truncation: solid 
- complete result, suming Eqs. (I22ap : dashed 
~ handbag contribution usually retained, 
Eq. (|17p : and dot-dashed - amendment 
described in connection with Eq. (|18l) . An 
interesting comparison is provided by the 
dotted curve: q 2 {x) = 30x^(l — To 

the eye, it is barely distinguishable from 
our complete result; and the mean value of 
the absolute relative difference between the 
curves is just 9%. 


the domain [0,1], and two convolution parameters z,z' G [—1,-|-1], so that the momentum 
integrals could be computed analytically. Inspecting the result, one can then determine Radon 
amplitudes for use in Eq. ([1]) that are consistent with those moments: the amplitudes vanish 
outside n = {(a,/3) : |a| + \j3\ < 1}, F{a,f3,t) is an even function of fd and G{a,fd,t) is odd. 
Eq. ([1]) then entails that H^{x,^,t) complies with the known constraints on polynomiality in 
vanishes outside x G [—^, 1] and is continuous at x = .^. The same recipe is to be applied 
to Eg. ([2nP in order to determine the correction’s Radon amplitude F^{a, (d,t). It worth to 
highlight, once again, that the results from Eqs. (11411 and (120^ . in the appropriate forward limit, 
exactly correspond with those from Eqs. (fT7P and (fTHp . respectively, shown by Eqs. (I22ap . 

We note now that when considering the comprehensive GPD defined by Eq. ([2]), one may 
write with complete generality: 


H^{x,0,-t) = i4^(x,0,0)AA(t)C^(x, t)F^(t), (23a) 

l=AA(t) J dx F[T^{x,0,0)Cn{x, t), (23b) 

so that all (x, t) correlations in H.,^ are expressed by C 7 r(x, t), which is necessarily non-unity in 
any physical system [12]. It is plain from Eq. ([1]) that only F(a,/3,t) contributes when = 0. 



Figure 5. Pion valence dressed-quark GPD, 
H^{x,0,—t), defined by Eos. (f23]l . ([2^ . 
obtained at the model scale, (h = 0.51 GeV; 



Figure 6. The same GPD evolved to 
C 2 = 2 GeV using leading-order DGLAP 
equations. 
















Summing the contributions from Eqs. (I14p . (120^ . the net result has the form 


F(q;,/3, -t) = 4>{a,P, tf[Fs{a,/3) + tV{a, P)(l){a, P, t)], (24a) 

4>{a, (3, t) = 1/[1 + (t/[4M2])(l - a + (3){1 - a - P)] , (24b) 

where the Fs component yields H^{x, 0, 0) = Qyix) in Ref. [T] and that with V is responsible for 
all violations of Eq. (fTOj) . 



Figure 7. Pion electromagnetic form 
factor obtained from F[^{x,0, —t), dehned 
by Eqs. ([23]), ([2^ . which is deliberately 
consistent with the result determined using 
Eqs. ([HI), (f^ and associated definitions. 
The data are described in Ref. |57] . The 
most favourable comparison is obtained with 
M = 0.40 GeV in Eqs. (l2T]l and the band 
shows results with M = 0.40 ± 0.05 GeV. 


Acting upon these observations, we define an ameliorated RL approximation to the pion’s 
valence dressed-quark GPD as the function obtained by: setting V = 0 in Eq. (I24a|) : and, for 
added simplicity, working with 4>{a, /3 = 0, f) whilst keeping the form of ^(^(t) computed directly 
from Eq. (1141) . Namely, via Eq. ([1]), our valence-quark GPD is given by Eq. ()23p with 

C(x,t) = l/[l + (t/[4M2])(l-x)^]2. (25) 


Our computed GPD is depicted in the upper panel of FigjSl Notably, the properties described 
in association with Eqs. (I10l) - (ll2p are evident, and this GPD naturally reproduces the pion 
valence dressed-quark distribution function obtained in Ref. [T]. 

The pion form factor associated with our GPD is drawn in Fig. [71 A ht to the result is 
provided by 

F^{t = M^z) = - 1 + 0.16Z - 

^ ^ 1-F 0.44z0.060^20.00033^3 ^ ^ 

At large-1 it behaves as 1/t^, whereas the correct power-law dependence is 1/t [58l [28l [29] . The 
power-law is wrong because Ea. (|21cp omits terms that have been described as representing the 
pion’s pseudovector components [59|, which are necessarily nonzero in a complete picture of the 
physical pion [231 [Hj. Notwithstanding that, it is valid and useful to compare the prediction 
with contemporary data and thereby determine a sensible value for our model’s dressed-quark 
mass-scale: the best comparison is obtained with M = 0.4 GeV. Notably, this scale is typical of 
the dressed-quark mass function in QGD [601 EH ES]- 

The IPD GPD derived from H^{x, 0, — t) in the upper panel of Fig. [5] is depicted in the upper 
panel of Fig. [HI The global maximum in this valence distribution is located at {x = 0.76, |6 _l| = 0) 
and, plainly, the probability to find a dressed-quark is strongly localised around this maximum. 
Naturally, for this valence dressed-quark distribution {d^\b_\_\ = 27rd|6_L| |&±|): 


f-idxf^ (f\b±\xql{x,\b±\) = 


1 

2 - 


(27) 














Figure 8. Pion’s valence dressed- 
quark GPD in impact parameter space, 
q^{x,\b±\](). Result obtained from 
HU in the top panel of 

Fig. [5] using Eq. (fT^ . 


Figure 9. Analogous result associated 
with H^{x, 0, — t; C 2 ) hr the bottom panel of 
Fig.O see Sect.5. [N.B. 1/M « 0.5fm, so 
b±M = 0.5 corresponds to b± ~ 0.25 fm and 
qlix,\b_L\]C)/M'^ = 1 means |6j_|; C) ~ 
4fm“^. 


5. Evolution of the GPD 

As explained elsewhere [T], our framework yields a valence-quark GPD that may be associated 
with an hadronic scale Ch = 0.51 GeV. It is worth outlining how the features of this distribution 
evolve to higher scales. Given that we have used ^ = 0, that is readily accomplished by using the 
DGLAP evolution equations to determine the evolved x-profile at each value of f. Our aim is to 
provide a qualitative illustration so, unlike Ref. pQ, we do not augment the valence distribution 
via the inclusion of gluon or sea-quark contributions. If desired, one could mask the impact of 
this omission by focusing on the behaviour of xH^{x,0,t) and xq^{x, |5_l|). 

Beginning with the valence dressed-quark GPD in the ufmer panel of Eig.El we employed 
leading-order evolution to obtain H^{x,0,—t;C2 = 2GeV)PJ The result is depicted in the 
lower panel of Eig.[5l Evidently, evolution, which adds glue and sea-quarks to the system by 
exposing this substructure within the dressed-quark, sharpens the peak associated with the 
global maximum at t = 0 and shifts its location toward x = 0. The maximum value at each 
t / 0 is also shifted toward x = 0; but outside a neighbourhood of t = 0 the profile in x is 
progressively flattened with increasing t. Notwithstanding this, at any finite C > C 2 , there will 
be a fi such that Vt > b a peak, albeit with much suppressed height, may be said to exist in 
the neighbourhood x ~ 1: ti increases with 

The last observation leads us to consider the conformal limit of QCD, which is recovered on 
r ~ 0, r = Aqcd/C- Within this domain, the valence dressed-quark GPD is i/^(x,0,0;r ~ 
0) = 6{x) [MIESIES]. [Eig. [5] highlights that this limit is reached slowly because evolution is 
logarithmic in QCD.] Eqs. ()23h then entail: 

HUx,0,-t;TC:^0) = 6{x)FUt). (28) 

This is a feature of our approximation to the pion’s valence dressed-quark GPD. It is not a 
property of the pion’s complete GPD, H-,^{x, 0, t; r ~ 0), because the valence GPD is a negligible 

^ Any significant differences generated by next-to-leading-order evolution are masked by a 25% increase in 
1 6.31 and hence are immaterial herein. 









piece of the full GPD on r ~ 0. That may be seen, e.g., by noting that valence-quarks carry 
none of the pion’s momentum within the conformal domain and hence it is invalid therein to 
represent Tn-(t) by an impulse (rainbow-ladder) approximation expressed through the triangle 
diagram of Eq. (fTTll . 

Having determined —t] C2), it is straightforward to obtain q^{x,\b±\-,(2) from 

Eq. ()13p . The result is depicted in the lower panel of Eig.[8l apparently, the maximum is shifted 
toward x = 0 and compressed in that direction, the peak height is diminished, and the width of 
the distribution in |6_i_| is increased. 

Each one of these evolution-induced changes may be intuitively understood by reasoning 
as follows. First consider a limiting case of an active parton with x ~ 1. This parton carries 
(almost) all the longitudinal momentum of the hadron. It therefore defines the CoTM and hence 
cannot be far removed from that centre. The distribution associated with an x ~ 1 parton must 
therefore be tightly localised around \b_i_\ =0. On the other hand, consider the case of an active 
parton with x reduced toward the location of the global maximum. The remaining partons 
within the hadron share in defining the CoTM and hence the active parton is not constrained to 
lie at |6_i_| =0. Plainly, as a parton’s value of x diminishes toward the favoured value, it plays 
less of a role in determining the CoTM and may therefore possess even larger values of |6j_|. 

In the current context, recall that evolution exposes the glue and sea-quark content of a 
dressed-quark: its identity comes to be shared amongst a host of partons, so that the probability 
of any one parton carrying x ~ I is much diminished. It follows that the global maximum in 
|6 j_ 1; C) must move toward x = 0 with increasing and, simultaneously, that the distribution 
is broadened in |6x| on the remaining domain of material support. 



Figure 10. Distribution of pion’s mean- 
square transverse extent, Eq. (I2^ : (black) 
solid curve (|6±(x; ^ 2 )^); and (blue) dashed 
curve - (|f>±(x; The (green) dotted 

curve is the result obtained at Qh when 
the correlation function C-,r{x,t) in Eq. (f2S|l 
is neglected: comparison with the dashed 
curve shows that the product Ansatz 
H^{x,0,—t) = gy(x)F 7 r(t) is generally 

unreliable. 


The latter effect is illustrated in Fig. llUl which depicts 

{\b±ix- C)P) = d^\b±\q{x, |6±1; C) |b±P ; (29) 

i.e., the x-distribution of the pion’s mean-square transverse extent: under evolution, the 
transverse extent narrows at large-x and broadens at small-x. A little consideration reveals 
that the measure of the curves in Fig. [10] is independent of the scale C because evolution is an 
operation that preserves the area under H{x,0,t) at each t. In fact, using Eqs. ([7|), (fTdIl . (I2^ . 
one finds 


/ I roo 

dx{\b_i{x;C)\‘^) = d\b±\\b±f d:^{\b±\), (30a) 

poo 

4(|6±|)=/ dAAJo(|6x|A)T;(A2), 

Jo 


(30b) 




and hence, with in Eq. (I26I1 . (|6 _lP) = (0.52 fm)^. For the reasons just explained, this is also 
the value obtained with Eq. (l28t) . [Note that En-(t) = l/(l + tr^/6), where is the pion’s electric 
charge radius, yields (|5±P) = (2/3)r^ = (0.55fm)^ (empirically [G^ = 0.672 ± 0.008fm).] It 

is natural that the mean-squared transverse extent of the dressed-quarks within a pion should be 
commensurate with the length-scale associated with light-quark confinement realised through a 
violation of reflection positivity [see, e.g.. Sect. 2.2 in Ref. mi 

6. Conclusion and prospects 

We described a calculation of the pion’s valence dressed-quark generalised parton distribution 
(GPD), within the context of a rainbow-ladder (RL) truncation of QCD’s Dyson- 

Schwinger equations. This framework is useful at an hadronic scale because it provides a 
description of hadrons via a dressed-quark basis, the accuracy of which in any given channel 
is knowable a priori. Our analysis focused here on the case of zero skewness (^ = 0) but we 
have also shown in Ref. [2] that the symmetry-preserving character of the RL truncation is 
all one needs to demonstrate a known relationship between = ±l,t) and the pion’s 

valence-quark parton distribution amplitude. 

Drawing analogy with the pion’s valence dressed-quark distribution function, we argued that 
the impulse-approximation used hitherto to dehne the pion’s valence GPD is generally invalid 
owing to omission of contributions from the gluons which bind dressed-quarks into the pion. 
We used a simple correction [Ea. ()20h ]. valid in the neighbourhood of ^ = 0, t = 0, in order to 
identify a practicable improvement to the approximation for Hi^{x,0,t). Expressing the result 
as the Radon transform of a single amplitude, we were able to isolate and remove those terms 
which produce unphysical behaviour, such as violations of positivity by the ^ = 0 GPD. The 
resulting, ameliorated Radon amplitude yields a form for H^{x,0,t) [Eqs. ([25]) . ([2^ ] which is 
consistent with signihcantly more known constraints than is the result produced by the impulse 
approximation alone. 

The results obtained in this way for H^{x,0,t), q^{x, |6 _l|) [Figs. El [5] provide a qualitatively 
sound picture of the dressed-quark structure of the pion at an hadronic scale. Using leading-order 
expressions, we evolved these distributions to a scale C = 2 GeV. All features of the resulting 
valence quark GPDs may be intuitively understood and hence the distributions should serve as 
an elementary but reasonable guide in the planning and interpretation of relevant experiments 
at existing or anticipated |68j facilities, which could plausibly involve deeply-virtual Gompton 
scattering on pions in a nucleon’s meson cloud. 

Notwithstanding the simplicity of the framework employed herein, a merit of the approach is 
its potential to compute features of hadron GPDs on the valence-quark domain and relate them 
directly to properties of QGD. This capacity has already been demonstrated in the simpler 
case of the pion’s valence parton distribution function [T]. One may begin to realise that 
potential by using more realistic forms for the dressed-propagators and -vertices that appear 
in the RL truncation analysis and, perhaps more importantly, uncovering the amendment to 
impulse approximation which is required in order to extend the validity of the RL truncation to 
the entire kinematic domain of ^ and t. 
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